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ABSTRACT 



G\ 



Context. Mon R2, at a distance of 830 pc, is the only ultracompact H n region (UC H n) where the photon-dominated region (PDR) between the 
i i ' ionized gas and the molecular cloud can be resolved with Herschel. Therefore, it is an excellent laboratory to study the chemistry in extreme PDRs 
<£ \ (G >10 5 in units of Habing field, n>10 6 cm" 3 ). 

OAims. Our ultimate goal is to probe the physical and chemical conditions in the PDR around the UC Hn Mon R2. 
Methods. HIFI observations of the abundant compounds "CO, C ls O, o-H 2 ls O, HCO + , CS, CH, and NH have been used to derive the physical 
(— i and chemical conditions in the PDR, in particular the water abundance. The modeling of the lines has been done with the Meudon PDR code and 
Qj, the non-local radiative transfer model described by Cernicharo et al. (2006). 

I ■ Results. The 13 CO, C ls O, o-H 2 18 0, HCO + and CS observations are well described assuming that the emission is coming from a dense 
O ' (ra=5xl0 s cm 4 , N(H2)>10 22 cirT 2 ) layer of molecular gas around the Hn region. Based on our o-H2 18 observations, we estimate an o-ITO 
J~i ' abundance of »2xl0 _8 . This is the average ortho-water abundance in the PDR. Additional H 2 18 and/or water lines are required to derive the 
C/3 , water abundance profile. A lower density envelope (n~10 5 cm" 3 , N(H2)=2-5xl0 22 cirT 2 ) is responsible for the absorption in the NH li— »02 line. 
C^ , The emission of the CH ground state triplet is coming from both regions with a complex and self-absorbed profile in the main component. The 
radiative transfer modeling shows that the 13 CO and HCO + line profiles are consistent with an expansion of the molecular gas with a velocity law, 
v e =0.5x(r/R o ,„)~' km s _I , although the expansion velocity is poorly constrained by the observations presented here. 

Conclusions. We determine an ortho-water abundance of »2xl0~ 8 in Mon R2. Because shocks are unimportant in this region and our estimate is 
JZi . based on H 2 18 observations that avoids opacity problems, this is probably the most accurate estimate of the water abundance in PDRs thus far. 

C^l ■ Key words. ISM: structure - ISM: kinematics and dynamics - ISM: molecules - HII regions - Submillimeter 

m ■ 

tT~ — ■ ; 1 . Introduction ated to IRS1 (Massi, Felli, & Simon 1985; Henning, Chini, 

( ^2 & Pfau 1992) and strong shocks are currently not at work in 

O Ultracompact (UC) Hn regions constitute one of the earliest this ion (Berne et al 2009) p revious moleC ular observa- 

r-H , phases in the formation of a massive star and are characterized tions (Giannak lou et aL 1997; Tafalla et aL 1997; choi et 

> ■ by extreme physical and chemical conditions (G„> 10 5 in units al 2000; RizzQ et al 20032005) showed that the uc Hn re _ 

._j , of Habing field and n> 10 cm" ). Their understanding is im- ion [& located inside a cayi and bound b a dense molecular 

^ ■ portant for distinguishing the different processes in the massive M The k of ^ molecular rid (herea fter, MP) is lo- 

H ; star formation process and because they can be used as a tern- cated at an offset (+10V1(r) relative to the ak of the ionized 

plate for other extreme photon-dominated regions (PDRs) such (here after, IF). The molecular hydrogen column density to- 

us the surface layers of circumstellar disks and/or the nuclei of ward ^ Mp ^ 2 _ 6xl0 22 cm -2 The detection of the reactiye 

starburst galaxies. The UC H n Mon R2 is the only one that can ions CQ+ and HQC+ showed a dense photon . dominate d region 

be resolved with Herschel „...-. , Q7fi , (PDR) surrounding the UC H n region (Rizzo et al. 2003, 2005). 

Mon R2 is a nearby (d=830 pc; Herbs & Racine 1976) Recent s -^ observations bed the thin mo i eC ular gas layer 

complex star forming region. It hosts a UC Hn region near its (n=4xl0 5 cm ^ N( h 2 )=1 X 10 21 cm" 2 ) with T,=574(±20) K in 

center, powered by the infrared source Mon R2 IRS1 (Wood between ^ ionized and ±e dense molecular traced b 

& Churchwell 1989). The molecular content of this region has ious millimeter obserV ations (Berne et al. 2009). All these 

been the subject of several observational studies. The huge co ents are schematically shown in Fig. 1. 

CO bipolar outflow (Meyers-Rice & Lada 1991), -15' long F J 6 
(=3.6 pc) is a relic of the formation of the B0V star associ- 

2. Observations 

* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- The observations were made with the HIFI instrument on- 
tant participation from NASA. board Herschel (Pilbratt et al. 2010, de Graauw et al. 2010) 
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Table 1. Summary of HIFI observations 



Fig. 2. HIFI spectra toward the MP in Mon R2. Note that the 
13 CO 5— >4 and HCO + 6— >5 lines appear close in velocity. This 
is an artifact of the DSB observations. The HCO + 6-5 line is 
observed in the LSB and the 13 CO 5->4 line, in the USB. The 
dashed line indicates v/ ir =10.5 km s . 

during the Priority Science Phase 2 in the frequency switch 
(FSW) observing mode with a reference position at the off- 
set (+10',0). Two receiver settings were observed, one in Band 
la with the WBS centered at 536.066GHz (LSB) and the 
other in Band 4a with the WBS centered at 971.800 GHz 
(LSB). These settings were observed toward both positions, 
IF [RA=06h07m46.2s, DEC=-06°23'08.3" (J2000)] and MP 
[RA=06h07m46.87s,DEC=-06°23'18.3" (J2000)]. In this letter 
we present the observations toward the MP because the study 
of the spatial distribution of the molecular tracers is postponed 
for a forthcoming paper. The data were reduced using HIPE 3.0 
pipeline. The adopted intensity scale was antenna temperature. 
In addition to Herschel data, we used the HCO + (l->0), HCO + 
(3->2), CS (2->l), I3 CO (l->0), 13 CO (2->l), C I8 (1^0) and 
C 18 Q (2->l) observed with the IRAM 30m telescope. 



10, o-H* 8 



3. Results 

We detected the CH 1^0, HCO + 6^5, CS 11 - ,.,,. 

lio^loi, 13 CO 5^4, C I8 5^4, and NH lj^0 2 lines in the 
MP (see Fig. 2). We stress that this is the first detection of the 
rarer water isotopologue o-H^O toward a spatially resolved 
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PDR. All the lines except NH 1 1 — >02 were detected in emis- 
sion. Gaussian fits are shown in Table 1 and Fig. 2. We fitted the 
three components of the CH 1 — »0 line assuming the same ex- 
citation temperature and estimated that the opacity of the main 
component is >3. In this case, the central velocity is not well 
determined because of the self-absorption and the flattened pro- 
files produced by the large opacities. The NH li— »02 line was 
tentatively detected (~3cr) in absorption. This line is composed 
of 10 components and we fitted all of them assuming the same 
excitation temperature. Our fit shows that the NH line is opti- 
cally thin. Because the individual components are not resolved 
with a linewidth of ~4 km/s, the individual line parameters are 
uncertain. This detection needs to be confirmed. 

3.1. Line profiles 

Different velocity components can be distinguished in this re- 
gion. The ambient cloud is centered at v/ 5r =10.5±l km s , A 
large scale molecular outflow is associated with IRS 1 (Meyers- 
Rice & Lada, C.J. 1991; Tafalla et al. 1997). The wings ob- 
served in the emission profiles of the HCO + 1— >0 and 3— >2 lines 
and in the 13 CO 2— »1 line (velocity ranges [0,6] km s -1 and 
[4,14] km s _1 ) are associated with the molecular outflow. The 
HCO + 6->5, 13 CO 5->4 and C 18 5->4 fines are centered at a 
velocity of ~10.4 km s _1 , In Fig. 3 we compare the HIFI lines 
with the low rotational lines of the same species observed with 
the IRAM 30m telescope. The angular resolution of the IRAM 
data was degraded to match those of Herschel. The line profiles 
of the low rotational lines of HCO + and 13 CO are self-absorbed 
at red-shifted velocities. For C ls O, there is a perfect match be- 
tween the profiles of the the J =5— >4 and J =2— >1 lines. The pro- 
files of the CS J=2— >1 and J=l 1— »10 lines also match perfectly. 
The HCO + 6^5, CS 11->10, 13 CO 5^4 and C I8 5->4 
lines have characteristic linewidths of 2-3 km s -1 . The largest 
linewidths observed in the NH 1 1 — >02 and CH lines (~ 5 km s -1 ) 
indicate that these lines are tracing a different, probably more 
diffuse component. As argued in Sect 4.1, the NH absorption is 
very likely caused by the cold and lower density envelope sur- 
rounding the UC H ii region. The CH 1 — »0 line is seen in ab- 
sorption and emission suggesting that CH is present in the low 
density envelope and the dense PDR. 

3.2. Molecular column densities 

In Table 2 we present the estimated molecular column densities. 
We used the rotational diagram technique to estimate the rotation 
temperature and the column density of the CO isotopologues and 
CS. This technique gives an average rotation temperature and to- 
tal column density providing that the emission is optically thin. 
For optically thick lines, we obtained a lower limit to the true 
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Table 2. Molecular column densities* 
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Fig. 3. Comparison of the profiles of the HIFI and IRAM 30m 
spectra toward the MP in Mon R2. The high angular resolution 
30m spectra were degraded to match the angular resolution of 
Herschel data (40"). The 13 CO l-»0 and C 18 l-»0 lines are 
not included because we observed them with a poor spectral res- 
olution (1.5 MHz~7 km s" 1 ). 
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* The rotational diagrams used the 30m and HIFI lines. 
" Assumed rotational temperature. 

* Assuming 12 C/ 13 C=50 and 16 O/ 18 O=500 
c From Rizzo et al. (2005) 

d Assuming an ortho-to-para ratio of 3. 



column density and rotation temperature (this is very likely the 
case for 13 CO). For CH and NH we assumed a reasonable value 
of the rotation temperature. Previous molecular observations at 
millimeter wavelengths showed the existence of gas with den- 
sities up to n=4xl0 6 crrT 3 and Tt^50 K in the molecular ridge 
(Choi et al. 2000; Rizzo et al. 2003, 2005). Berne et al. (2009) 
derived a density, n=4xl0 5 crrT 3 and a gas kinetic temperature 
of ~570 K for a thin gas layer around the H n region on basis of 
the H2 rotational lines. These two n-T pairs of values can pro- 
duce the rotation temperatures derived from our observations. 
However, the gas layer at 570 K cannot account for the observed 
line intensities. For instance, assuming a standard C 18 abun- 
dance of 2xl0~ 7 , this hot layer will produce a C ls O 5— A line 
intensity of T/,~0.09 K, i.e. 2% of the observed value. Therefore, 
the emission of the observed lines is mainly coming from the 
dense (n=4xl0 6 cm 4 ) molecular ridge. 

We detected the o-H^O li ->loi line in the MP. Using 
the non-local radiative transfer code of Cernicharo et al. (2006) 
and assuming the densities and temperatures prevailing in the 
dense molecular ridge (Tt=50 K, n(H2)=4xl0 6 crrT 3 ), we ob- 
tained an excitation temperature for the ground state transition 
of o-H^O of ~8.5 K. With this low excitation temperature, 
we need N(o-H2 8 0)= 2.7xl0 12 cirT 2 to fit our observations. 
Assuming a ls O/ I6 ratio of ~500, this implies an ortho-water 
abundance of ~2xl0~ 8 . The excitation temperature is not very 
sensitive to the gas kinetic temperature. To assume a gas kinetic 
temperature as high as 500 K would increase the excitation tem- 
perature by a factor of 2, and decrease the estimated o-H^O 
column density by a factor of 10. To assume a lower hydrogen 
density would be more critical, because the excitation tempera- 
ture would drop to very low values (<6 K), and the line would 
become very weak and optically thick which would prevent any 
good estimate of the o-H^O column density. The water abun- 
dance estimated in Mon R2 is similar to that obtained toward the 
Orion Bar by Olofsson et al. (2003) using ODIN observations. 
In that case, the spatial resolution of the ODIN observations did 
not allow the authors to resolve the dense PDR. Moreover, their 
estimate was based on observations of the ground state transi- 
tion of the main water isotopologue, which is optically thick. 
The agreement between the two measurements could therefore 
be fortuitous. Recent SPIRE observations of the Orion Bar have 
provided an upper limit to the water abundance of a few 10~ 7 
(Habartetal. 2010). 

To interpret the observed NH 1 1 -O2 line absorption, a knowl- 
edge of the continuum level is needed. Unfortunately, the line 
was observed in FSW mode, which removes the continuum. 



Hence, the continuum level had to be estimated from previous 
measurements. In particular, Dotson et al. (2010) measured an 
intensity of 280 Jy/beam using the CSO telescope (beam=20") 
at 350 /im (=857 GHz). Assuming a spectral index/? of ~3, we 
estimate a continuum flux of 411 Jy at 974 GHz. This value 
corresponds to a brightness temperature of 0.8 K in the HIFI 
beam. Taking into account the uncertainty in /3 and the differ- 
ent beams of Herschel and CSO, we estimate that the accuracy 
of the continuum intensity at 974 GHz is about a factor of 2. 
Assuming T ex - 10K, and an intrinsic linewidth of 4 km s -1 , the 
line opacity inferred from the observed absorption feature im- 
plies a NH column density of (l-5)xl0 13 cm 2 (a factor ~1.5 



higher if T ev =20K). Taking ~5xl0 cm as an upper limit 
for N(H2) (the NH absorption likely arises in a external layer 
of more diffuse gas), the NH abundance would be greater than 
(1.0-0.2)xl0~ 9 . These abundances are comparable to the NH 
abundance first inferred by ISO toward Sgr B2 (Cernicharo et 
al. 2000, Goicoechea et al. 2004). 



4. Discussion 

4.1. Chemical model 

We explored the possibility of explaining the molecular abun- 
dances observed in Mon R2 in terms of PDR chemistry. To do 
this, we used the updated version of the Meudon PDR code 
(Le Petit et al. 2006, Goicoechea & Le Bourlot 2007). As in- 
put parameters we used a plane-parallel slab with a thickness 
of 10 mag and n=4xl0 6 cirT 3 , which is illuminated from the 
left side with a field of Go=5xl0 5 and from the right side with 
the standard interstellar UV field (Go=T). Elemental abundances 
were: He (0.1), O (3.19xl0~ 4 ), N (7.5xl0 5 ),C (1.32xl0~ 4 ), S 
(1.86xl0~ 5 ), Fe (1.5xl0~ 8 ). The model does not include iso- 
topic fractionation, we estimated the abundances of the iso- 
topologues assuming 12 CO/ 13 CO=50 and 16 O/ 18 O=500. The ob- 
tained column densities are shown in Table 2 and the abundance 
profiles are plotted in Fig. 4. We got an excellent agreement be- 
tween observations and model predictions for the CO isotopo- 
logues, HCO + , H2O, and CS. The model successfully predicts 
the observed average water abundance. The water abundance is 
far from uniform across the PDR, with a peak of ~10~ 6 at an 
extinction of ~ 1 mag and a minimum value <10~ 10 (see Fig. 4). 
The observation of high excitation lines of water will allow us 
to derive the water abundance profile and further constrain the 
chemical modeling. The model falls short by more than one or- 
der of magnitude however to predict the NH column density. The 
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NH abundance is <10~ 9 all across the PDR (see Fig. 4). As ar- 
gued below, this is very likely due to the main contribution to the 
observed NH line coming from the low density envelope that is 
protected from the UV radiation. We ran the code for a low den- 
sity (nx4xl0 4 cm 4 ) plane-parallel layer of 50 mag illuminated 
with Go=l and obtain a NH column density of ~3xl0 13 crrT 2 in 
better agreement with the observations. This low density enve- 
lope would also contribute to the CH ground state line since the 
CH column density should be ~ a few 10 14 cm 2 . For this reason, 
we have a very complex CH profile, with the line optically thick 
and self-absorption in the main component. However, the con- 
tribution of the low density envelope to the other observed HIFI 
lines would be negligible because of the low density and low gas 
kinetic temperature. For instance, the intensity of the 13 CO 5— >4 
line would be ~0. 1 K. 



4.2. Kinematics of the region 

The systematic change in the line profiles of the rotational lines 
of HCO + and I3 CO can be understood in terms of the velocity 
structure of the molecular core. Although the complete model- 
ing of the source is beyond the scope of this paper, (most of 
the HIFI data are still to come), we made a preliminary model 
to gain insights into the kinematics of this source. The UC H n 
region is expected to expand because of the different pressure 
of the ionized and molecular gas (Jaffe et al. 2003; Rizzo et al. 
2005). We modeled the HCO + and 13 CO lines using the non- 
local radiative transfer code of Cernicharo et al. (2006). In the 
model, we adopted the physical structure derived in Sect. 4.1 (a 
spherical envelope with inner radius R,„=0.08 pc, an innermost 
layer with Tt=50 K, n=4xl0 6 cirT 3 and a radial thickness of 
0.0006 pc, and an external envelope with the temperature de- 
creasing as TxR ° 5 , a constant density of n=lxl0 5 cirT 3 and 
a thickness of 0.16 pc). For the dust temperature and opacity, 
we adopted the values derived by Thronson et al. (1980). As a 
first approximation we assumed constant 13 CO and HCO + abun- 
dances of 2xl0~ 6 and 10~ 9 respectively and only varied the ve- 
locity law. We were able to reproduce the line intensities and the 
trend observed in the line profiles assuming an expansion veloc- 
ity law v e = 0.5x(R/R o „,) _1 km s _I (Ro Uf =0.24 pc) and a turbu- 
lent velocity of 2 km s (see Fig. 5). Because of the hole in the 
molecular core, the maximum expansion velocity is 1.5 km s -1 , 
which is similar to the turbulent velocity. For this reason, the 
back and front parts of the envelope are radiatively coupled pro- 
ducing the self-absorbed profiles at red-shifted velocities. The 
shapes of the line profiles are dominated by the turbulent ve- 
locity with little effect of the small expansion velocity, which is 
poorly constrained. 



5. Conclusions 

We present the first HIFI observations toward the UC Hn 
region Mon R2. Detections of the 13 CO 5^4, C 18 5->4, 
o-H 2 18 1io->loi,HCO + 6^5,CS 1 1— >10, NH l r 2 , and CH 
1— >0 lines are reported. The emission/absorption of all these 
molecules is well explained assuming the the molecular core is 
composed of a dense (n=5xl0 6 cm 3 ) PDR layer of gas sur- 
rounded by a lower density UV protected envelope. The mod- 
eling of the 13 CO and HCO + line profiles is consistent with the 
molecular gas being expanding with an expansion velocity law, 
v e =0.5x(r/R„„,)~' km s -1 . Based on our o-H 2 18 lio->loi ob- 
servations, we estimate an ortho-water abundance of ~2xl0~ 8 . 
Because shocks are unimportant in this region and our estimate 
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Fig. 5. Predicted line profiles for the 13 CO and HCO + rotational 
lines using the model described in Sect 4.2. 



is based on the rarer isotopologue observations that avoids opac- 
ity problems, this water abundance estimate is probably the most 
accurate in PDRs thus far. 
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Fig. 1. Overview of the PDR associated with Mon R2. Gray 
scales represent the CS 7— »6 line emission map from Choi et al. 
(2000). The blue contours represent the intensity of the [Nen] 
line (0.02-0.1 erg s cirT 2 sr~' in linear steps), and the green 
contours the intensity of the H2 S(3) rotational line at 9.7 fim 
(1.5^4.5 xlO 4 erg s~ : cirT 2 sr _1 in linear steps). The circles in- 
dicate the Herschel beam in band la (40") and band 4a (21") 
centered on the molecular peak (MP). 
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Fig. 4. Results of the chemical modeling of the dense PDR sur- 
rounding the UCHII region Mon R2 



